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Abstract 

Events from the e"'"e~ — > Z7 process with hard initial-state radiation collected 
with the L3 detector at centre-of-mass energies between 183 GeV and 209 GeV are 
used to measure the mass of the Z boson. Decays of the Z boson into hadrons or 
muon pairs are considered and the Z mass is determined to be 91.272±0.032 (stat.)± 
0.033 (syst.) GeV, in agreement with the value measured at the Z resonance. Alter- 
natively, assuming this measured value of the Z mass, the method determines the 
LEP centre-of-mass energy, found to be 175 ±68 (stat.) ±68 (syst.) MeV lower than 
the nominal value. 



Submitted to Phys. Lett. B 



1 Introduction 



At the LEP collider, operating at centre-of-mass energies, a/s, above the Z peak, the process 
e"'"e~ — > ff7 frequently occurs with hard initial-state radiation (ISR). Due to the strong reso- 
nance in the reaction e"'"e~ — > Z — ff , the energy of the ISR photon is such that the mass of 
the fermion pair recoiling against the photon is often close to the mass of the Z boson, mz- The 
strong forward peaking of the ISR process yields photons which usually remain undetected in 
the beampipe but in a small fraction of events they are detected. 

Data recorded with the L3 detector [1] at ^/s = 183 — 209 GeV with a total integrated 
luminosity of 685 pb~^, detailed in Table [TJ are used to extract the mass of the Z boson, 
using Z decays into quark or muon pairs. The mass of the hadronic system after applying a 
kinematic fit, rUes, is directly related to mz- In the case of muon pairs the event kinematics is 
fully determined by the measurement of the muon scattering angles, making use of the excellent 
polar angle resolution of the muon spectrometer. From the measured or reconstructed photon 
energy, E^, the effective centre-of-mass energy of the muon pair, y/s', is calculated as 

An unbinned likelihood fit is applied to the distribution of nies or to extract the Z mass. 
The value is compared with the precision measurement of L3, derived from the Z-lineshape scan 
at centre-of-mass energies around the Z pole [2] . This comparison serves as a cross check of the 
W-mass measurement which uses similar techniques. Assuming the Z mass measured around 
the Z pole, the method provides a measurement of the LEP centre-of-mass energy which is 
compared to the one determined by the LEP energy working group [3]. 

2 Monte Carlo simulation 

The Standard Model predictions for the different final states are determined with the following 
Monte Carlo programs: KK2F [4] for e+e~ —>■ qq(7), e+e" — *• /i+/x~(7) and e"'"e~ T~^T~{'y), 
BHWIDE [5] for e+e" ^ e+e-(7), PHOJET [6] for e+e" ^ e+e'hadrons, DIAG36 [7] for 
e"'"e~ — ^ e~^e~fi^fi~, KORALW [8] for e"'"e~ W+W^, with the exception of the ez/qq' final 
state, which is modelled by EXCALIBUR [9], and PYTHIA [10] for e+e" ZZ and e+e" ^ 
Ze"'"e~. The response of the L3 detector is modelled with the GEANT [11] detector simulation 
program which includes the effects of energy loss, multiple scattering and showering in the 
detector material. The GEISHA program [12] is used to simulate hadronic interactions in the 
detector. Time-dependent detector inefficiencies are taken into account in the simulation. 

3 Selection of hadronic events 

To remove purely leptonic final states the event must have more than 12 calorimetric clusters. 
In addition, the transverse energy imbalance must be less than 25% of the visible energy and 
the sum of the cluster energies must be greater than 30% of a/s, as shown in Figure 

Hadronic final states from two-photon collisions are typically boosted along the beam di- 
rection. To reject these events we require the longitudinal energy imbalance not to exceed 80% 
of the visible energy. Additionally, the sum of the transverse energies of the hadronic clusters 
has to be greater than 15% of ^/s. Both cuts effectively remove the two-photon background. 
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Four-jet final states from W- and Z-pair production have a total hadronic energy close to 
^/s. The clusters in these events are distributed more spherically than in the case of two-jet 
events from Z7 production. For events with a sum of the hadronic energy greater than 70% of 
^/s, the particles are boosted to the rest frame of the hadronic system. A cut on the thrust, T, 
of the boosted event, T > 0.85, removes four-jet events as shown in Figure^. 



4 Selection of muon-pair events 

Muons are reconstructed from tracks in the muon chambers in a fiducial volume of | cos6\ < 0.9. 
If only one muon track is present in the event, the signature of an additional minimum- ionizing 
particle in the inner detectors is required. Hadronic events are suppressed by requiring less than 
15 calorimetric clusters. To reject cosmic-ray background, the time measured by one scintillator 
matched in azimuthal angle to a muon track has to be consistent with that of beam crossing 
within 5 ns. In addition, the distance of closest approach of at least one of the muon tracks to 
the interaction point in the plane perpendicular to the beam must be less than 1 mm. 

The two muons are usually almost back-to-back in the plane perpendicular to the emitted 
photon. Therefore, the angle between the muon directions in that plane, A0, is required to be 
greater than 175°, as shown in Figure The measured momentum, p^, of the muon which 
has the largest polar angle, 9i, has to be greater than 60% of the corresponding expected 
momentum, p'j^^, calculated according to 

exp _ /- sin 6*2 , s 

~ sin 01 + sin 02+ I sin (01 + 02)1' ^ 

where 02 is the polar angle of the other muon. Figure^ shows the distribution of Pfj,/pJ^^- 



5 Mass reconstruction 

In each selected event we search for an isolated photon. A photon candidate is a calorimetric 
cluster with an energy greater than 30% of ^/s. It has to be isolated, with at most one additional 
cluster within 15° of the candidate direction. In addition, the angle between the photon and 
the nearest track must be greater than 4.6°. 

Hadronic events are forced into two jets using the DURHAM [13] algorithm, excluding 
photon candidates, if any. The pion mass is assigned to energy depositions matched with a 
charged track while clusters without a matching track are treated as massless. A kinematic fit 
is applied to the event, where the measured quantities are varied within their resolution and 
four-momentum conservation is imposed. If no high-energy photon is observed in the detector, 
a single photon is assumed to have escaped undetected along the beam direction. This improves 
the resolution of the two-jet effective mass by a factor of about 3. The distribution of the fitted 
mass, rries, is shown in Figure |21 

In muon-pair events, the photon energy is calculated from the reconstructed muon momenta 
using three-particle kinematics as 

E = |sin(0i + 02)| 

^ ^ sin0i + sin02+ |sin(0i + 02)1 ■ 

The 0j are the angles between the muons and the photon direction. If no photon is detected 
they are the angles between the muons and the beam direction. The effective centre-of-mass 
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energy of the muon pair, v^, is calculated using Formula [T] and its distribution is shown in 
Figure El 



6 Fit method 

Events with hard ISR are selected by requiring 

70 GeV < meff < 110 GeV 

for hadronic events and 



80 GeV < Vs' < 100 GeV 

for muon-pair events. Data samples of 34081 and 799 events are selected for hadrons and muon 
pairs, respectively, with background levels of 3.2% and 6.2%. 

To extract the Z mass from this data sample, an unbinned likelihood fit is applied to the 
measured differential cross section as a function of nies or The likelihood is defined as the 
product of the normalised differential cross sections of each event after selection cuts: 

where m|* is the Z mass varied during the fit and ^ represents rrics for the hadronic events and 
y/s' for the muon-pair events. The total and differential accepted cross sections of the signal 
are denoted as a and da/d^, while ctbg and duBc/d^ are the total and differential cross sections 
of the background. 

The box method [14] is used to obtain the accepted differential cross section. This method 
takes into account both detector resolution and selection efficiency effects. The accepted dif- 
ferential cross section is determined by averaging Monte Carlo events inside a ^ bin centered 
around each data event. The number of signal Monte Carlo events is scaled such that it agrees 
with the measured total cross section, while the number of background events is normalised 
to the integrated luminosity. This ansatz assigns fluctuations in the number of data events to 
the signal. Stable results are obtained with bin sizes chosen such that 500 signal Monte Carlo 
events are contained in each meg bin and 250 signal events in each y/s' bin. The size of the bins 
ranges from about 20 MeV in the rries peak to about 6 GeV in the rries tails and from 50 MeV 
in the y/s' peak to 3 GeV in the y/s' tails. 

To simulate the effect of different values of the Z-boson mass with a finite number of Monte 
Carlo events, the signal Monte Carlo events are reweighted using the weights 

^r,<') . (5) 

where yfs' is the generated effective centre-of-mass energy of each event and is the value 
of the Z mass that is used in the Monte Carlo generation. The calculation of the differential 
cross section, da/dx/s', up to 0{a'^) of Reference 15 is used. 

To confirm the linearity of the fitting method as well as to test for any possible bias, five 
Monte Carlo samples, each with about ten times the statistics of the data, are generated with 
different Z masses between 90 GeV and 92 GeV. The correct values for the mass of the Z boson 
are found by the fitting method within the statistical precision of the test. 
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7 Systematic uncertainties 



The systematic uncertainties for the Z mass determination are hsted in Table |21 and are detailed 
below. 

To estimate the uncertainty due to hadronisation, three different Monte Carlo samples 
are generated which have the same underlying two-fermion events, but different hadronisation 
schemes, modelled by the programs ARIADNE [16], HERWIG [17] and PYTHIA. As the jets 
are formed from clusters that are either massless or assigned the pion mass, a different kaon 
or baryon content in data and Monte Carlo would lead to a bias in the extracted Z mass. 
For each of the three models the Monte Carlo events were reweighted in order to reproduce 
the mean kaon and proton multiplicity measured on the Z peak [18]. The mass of the jets 
plays a role in the kinematic fit and the differences of the jet masses produced by the different 
hadronisation schemes are also considered. A reweighting method is applied to reproduce the 
jet masses measured in data. The root-mean-square of the various Z masses obtained with 
the three Monte Carlo programs after the different reweighting procedures is assigned as the 
hadronisation systematic uncertainty. 

The ratio of the measured jet energy to the energy calculated from the jet angles analogously 
to Formula El is shown in Figure EK. The jet energy scale is confirmed at the percent level. The 
value of this ratio as a function of cos 6 allows for the recalibration of the energy measurement 
of hadronic clusters. Half of the difference to the mass obtained without recalibration is taken 
as systematic uncertainty. In addition, we scale the energy of the individual energy depositions 
by ± 0.5%. The difference to the value without scaling is taken as systematic uncertainty. The 
total uncertainty from energy calibration is the sum in quadrature of the two contributions. 

The uncertainty due to the measurement of the jet polar angle is estimated by repeating the 
analysis with jets formed by calorimetric clusters that have a corresponding charged track and 
jets formed by the same clusters but using the polar angles from the correlated charged track 
instead of the angle measured in the calorimeter. The accuracy of the measurement of the muon 
polar angle is tested by comparing the default ^fs' spectrum to the spectrum obtained using 
the angles of calorimetric clusters associated with the muons. The differences of meff or \fs' are 
calculated on an event-by-event basis and shown in Figures lllo and lit, respectively. We assign 
half the difference of the average shifts of data and Monte Carlo as systematic uncertainties. 

Additional clusters around detected photons, wrongly assigned to one of the jets, would 
affect the reconstructed Z mass. Removing all clusters in a 10° cone around the photon direction 
results in a negligible mass shift. In a similar way, calorimetric clusters from random noise, 
equally distributed in the calorimeters, would affect the reconstructed mass if they are not 
described by the Monte Carlo. The angular distribution of calorimetric energy relative to the 
corresponding jet axis is plotted in Figure EH- Good agreement between data and Monte Carlo 
is seen. Removing all clusters outside a cone of 60° half-opening angle around the jet axis yields 
a negligible mass shift. 

For both hadronic events and muon-pair events the uncertainty in the background level is 
evaluated by scaling the total cross section of the background Monte Carlo samples by ±5%. 

The theoretical uncertainty on the ISR modelling is tested by comparing results obtained 
with different ISR modelling schemes to different orders of a. This is done by using event 
weights given by the Monte Carlo generator KK2F as described in Reference 4. The differences 
between the results from 0{a^) calculations with Coherent Exclusive Exponentiation with and 
without ISR/FSR interference is assigned as systematic uncertainty. 

The stability of the box method is tested by changing the box sizes used in the fit. A small 
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shift is observed and quoted as systematic uncertainty. To test the resuhs obtained with the 
maximum hkehhood fit, we also apply a fit of the Monte Carlo histogram to the measured 
meff distribution by reweighting individual Monte Carlo events. As a second test, we apply 
the maximum likelihood fit to the ^J~s' spectrum of the hadronic event sample, calculated using 
Formula El modified to include the effect of non zero jet masses. The results from both cross 
checks show no significant deviation from the default method. The uncertainty from limited 
Monte Carlo statistics is also considered. 

The LEP beam energy has an uncertainty of 10 to 20 MeV [3] depending on ^/s. The 
relative error on the Z-boson mass is the same as that on the beam energy. 

The systematic uncertainties are treated as fully correlated between the different energy 
points, except the uncertainty resulting from limited Monte Carlo statistics, which is treated 
as uncorrelated, and the uncertainties on the beam energy, where the correlation matrix from 
Reference [3] is used. A total systematic uncertainty of 39 MeV for the determination of the 
Z-mass from the hadronic channel and 21 MeV from the muon channel is found. 

8 Results 

The results for different values of \fs are shown in Table El They are combined taking the 
systematic uncertainties into account. The combination yields 

mf- = 91.271 ± 0.031 (stat.) ± 0.039 (syst.) GeV 

for hadronic events and 

mf^'' = 91.276 ± 0.105 (stat.) ± 0.021 (syst.) GeV 

for muon-pair events. The results of the fit are shown in Figures |21 and El Averaging the results 
obtained from the hadronic and muon pair samples, including all correlations, yields 

^meas ^ 272 ± 0.032 (stat.) ± 0.033 (syst.) GeV. 

For the combination of hadron and muon-pair events, the uncertainty due to ISR modelling 
is treated as fully correlated. For the beam energy uncertainties the correlation matrix of 
Reference 3 is used. This value is in good agreement with the precision mass measurement, 
mz = 91.1898 ± 0.0031 GeV [2]. 

The measurement can also be interpreted as a determination of the LEP centre-of-mass 
energy, ■^/s"^^'^. A difference between the measured mass, m^'^^, and the precision mass, mz, 
can be attributed to a deviation from the nominal value, ^/s■. 

Av/i = y^"^^^ - = -V~s^ ^. (6) 

mz 

The value obtained from the observed Z mass, 

Av^ = -0.175 ± 0.068 (stat.) ± 0.068 (sys.) GeV, 

is consistent with no shift. 

In conclusion, the Z mass measured in radiative events is in agreement with the determina- 
tion at the Z pole, validating the method used for the measurement of the mass of the W boson. 
Interpreted as a determination of the centre-of-mass energy it agrees with the measurements 
by the LEP energy working group [3]. 
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Year 


(yi) [GeV] 


Luminosity [pb ^ 




1997 


182.7 


55.5 


1998 


188.6 


176.8 


1999 


191.6 


29.8 




1999 


195.5 


84.1 




1999 


199.5 


83.3 




1999 


201.8 


37.2 




2000 


204.8 


79.0 




2000 


206.5 


130.8 




2000 


208.0 


8.3 




Total 




684.8 



Table 1: Centre-of-mass energies and corresponding integrated luminosities. 





Uncertainty [MeV] 


Uncertainty [MeV] 


Source 


(hadron channel) 


(muon channel) 


Hadronisation 


22 




Energy calibration 


16 




Angular measurement 


24 


11 


Background 


3 


<1 


Initial- state radiation 


4 


11 


Box Size 


2 


2 


Monte Carlo statistics 


8 


9 


LEP energy 


11 


11 


Total 


39 


21 



Table 2: Systematic uncertainties on the Z boson mass. 





(hadron chanuel) 


(niuou channel) 


182.7 
188.6 
191.6 
195.5 
199.5 
201.8 
204.8 
206.5 
208.0 


91.286 ±0.095 ±0.046 
91.290 ±0.057 ±0.042 
91.402 ±0.143 ±0.046 
91.467 ±0.089 ±0.046 
91.144 ±0.094 ±0.046 
91.073 ±0.142 ±0.046 
91.369 ±0.101 ±0.042 
91.107 ±0.081 ±0.042 
91.329 ±0.331 ±0.042 


91.057 ±0.317 ±0.029 
91.224 ±0.189 ±0.029 
92.065 ±0.635 ±0.029 
91.219 ±0.332 ±0.029 
91.183 ±0.422 ±0.029 
91.464 ±0.402 ±0.029 
91.358 ±0.260 ±0.032 
91.439 ±0.273 ±0.033 
90.439 ±0.667 ±0.033 


Combined 


91.271 ± 0.031 ± 0.039 


91.276 ±0.105 ±0.021 



Table 3: Results of the maximum hkelihood fit for each ^/s value. The first uncertainty is 
statistical, the second systematic. 
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Figure 1: Variables used for the event selections: a) visible energy, normalised to y/s, b) 
thrust in the centre-of-mass frame of the jets, c) angle between the muons in the 
plane perpendicular to the photon and d) measured muon momentum, normalised to 
the expected momentum. 



12 



3000 



> 2000 
CD 



if) 

miooo 







70 



Data 



qq(y) 



A 



L3 



Background 





80 



90 



n — ■ — T — ■ — j — ■ — I — — r — ' I ' ' 

100 110 



^ett [GeV] 



Figure 2: Distribution of the effective mass of the two-jet system after kinematic fit. The sohd 
hne shows the result of the fit which determines the Z boson mass. 
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Figure 3: Distribution of the effective centre-of-mass energy for muon-pair events. The sohd 
hne shows the result of the fit which determines the Z boson mass. 
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ure 4: Variables used for systematic studies: a) ratio of measured jet energy and jet energy 
calculated from jet angles, b) difference, per event, of the effective mass as obtained 
using cluster or track angles, c) difference, per event, of the effective centre-of-mass 
energy, obtained by using polar angles from the direction of the muons or the cor- 
responding calorimetric clusters, d) relative energy flow in a jet as a function of the 
angular distance from jet axis. 
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